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Abstract 

We show the traditional rocket problem, where the ejecta velocity is assumed 
constant, can be reduced to an integral quadrature of which the completely 
non-relativistic equation of Tsiolkovsky, as well as the fully relativistic equa¬ 
tion derived by Ackeret, are limiting cases. By expanding this quadrature 
in series, it is shown explicitly how relativistic corrections to the mass ra¬ 
tio equation as the rocket transitions from the Newtonian to the relativistic 
regime can be represented as products of exponential functions of the rocket 
velocity, ejecta velocity, and the speed of light. We hud that even low or¬ 
der correction products approximate the traditional relativistic equation to 
a high accuracy in flight regimes up to 0.5c while retaining a clear distinc¬ 
tion between the non-relativistic base-case and relativistic corrections. We 
furthermore use the results developed to consider the case where the rocket 
is not moving relativistically but the ejecta stream is, and where the ejecta 
stream is massless. 

Keywords: special relativity, relativistic rocket, Tsiolkovsky equation, 
Ackeret equation 


1. Introduction 

A fundamental result in the theory of rocket propulsion is the relation¬ 
ship between the velocity increment a rocket gains and the amount of reaction 
mass it expels to gain this increment. The canonical form of this equation. 
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which is credited to Tsiolkovsky, is derived under the assumption that the 
rocket is moving slowly compared to the speed of light throughout its flight 
([3]), and thus its momentum is given by the Newtonian formula p = rriR'v. 
Another relationship between rocket velocity and expelled reaction mass gov¬ 
erns the flight of rockets operating in the relativistic regime. The relativistic 
rocket equation was first derived by Ackeret ([2]) in 1946, and subsequently 
expounded upon by Krause m), Bade ([S]), and Forward mM), the last of 
whom presented a “transparent derivation” of the relativistic rocket equation 
based on simple equations in special relativity. Other research has been done 
in relativistic rocketry, such as simultaneity calculations based on the prop¬ 
agation of light signals to and from a constantly accelerating rocket m), 
calculations of the maximum attainable velocity for a relativistic rocket ([Sj) 
and how an Igp like parameter could be constructed for a rocket using exotic 
propellants whose products consist of both massive and massless particles 
( 0 , 0 ). 


One issue which has not yet been investigated is how a rocket will perform 
in a weakly relativistic regime, and how quasi post-Newtonian corrections can 
be derived for the Tsiolkovsky equation. In fact, while it has been noted by 
Forward ([T]) that the Ackeret rocket equation converges to the Tsiolkovsky 
equation if evaluated numerically, the explicit proof of this fact does not ap¬ 
pear to have been previously published. We shall show that starting from 
the same simple equations as Forward it is possible to reduce the rocket 
equation to a general quadrature, without making any assumptions about 
its flight regime. Moreover, this quadrature is shown to reproduce the equa¬ 
tions of both Tsiolkovsky and Ackeret when evaluated in the correct limits. 
Finally it is shown that by expanding this quadrature in series it is possible to 
obtain a product of quasi post-Newtonian corrections which converges to the 
Ackeret equation in the limit as n —)■ cx). Section 2 presents the derivation of 
the main result of this work. Equation (15), while Section 3 shows how the 
Tsiolkovsky and Ackeret equations as well as the post-Newtonian terms can 
be all derived from this result. Section 4 uses the results we develop to ana¬ 
lyze two seemingly novel scenarios, which are found to yield results that are 
surprisingly consistent with a naive application of the Ackeret’s relativistic 
equation. 
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2. The Rocket Problem and its Reduction to Quadrature 

In this section we discuss some basic relativistic principles necessary for 
the derivation and show how the quadrature is derived. 

2.1. Assumptions and Definitions 

Just as energy and momentum are conserved separately in Newtonian 
mechanics, energy-momentum is conserved as a 4-vector in special relativity. 
We dehne the energy-momentum of a system in (l-l-l)-dimensions as 

P/. : (Po,Pi) = (^/c,p). (1) 

Similarly, we may dehne a (l-l-l)-dimensional 4-velocity which is given by 

: (mo,Mi) = (7C,7'i'), (2) 

where 7 is the Lorentz factor of the frame which is dehned to be 



One may note that in this form the familiar Newtonian relations still hold, 
viz. 

Pf, = mUf,, (4) 

which the reader may verify. 

In order to perform any calculation in relativistic mechanics we must 
hrst specify a frame for the observer, which will be considered to be at rest. 
The most convenient frame for a rocket calculation will be the Center-of- 
Momentum (COM) frame, which is dehned as the unique frame for which 

= 0 - ( 5 ) 

j 

It can be shown in Newtonian mechanics that the center of mass of any closed 
system possesses a constant velocity. A similar conservation principle holds 
in special relativity and the COM frame may be thought of as the unique 
inertial frame which is comoving with the center of mass and whose origin 
coincides with the position of the center of mass. 
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Center-of-Momentum Frame 


Figure 1: Schematic of the energy-momentum balence. 


2.2. The Rocket Problem 

In the rocket problem, a body which is initially moving with velocity Vo 
(which is zero in the COM frame) ejects an inhnitesimal amount of mass 
dm which propels the body forward in accordance with the conservation of 
momentum. In the relativistic picture this is equivalent to the conservation 
of energy-momentum as has been dehned in Section 2.1. In order to form 
the energy-momentum balance we consider the rocket at two inhnitesimally 
separated instants: 

(i) : The rocket possesses zero velocity with respect to the COM frame and 
has rest mass m + dm. We shall consider dm to be the rest mass of the 
propellants so that after they are expelled the rocket possesses a rest mass 

m. 

(ii) : The rocket expels the mass dm and gains a velocity increment dw. This 
separates the center of mass of the rocket from the center of momentum of 
the system and thus gives the rocket a nontrivial velocity with respect to the 
COM frame. 
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The process is further illustrated in Figure 1. As stated before, we may write 
energy-momentum conservation in the COM frame as 

d = 0 . ( 6 ) 

For the rocket system the energy-momentum consists of two components, the 
rocket component and the ejecta component p®. Substituting this into 
Equation ([^ we hnd 

<ipj = -dPl- (7) 

Because the ejecta is presumed to be at rest to begin with, the differential 
momentum for the ejecta is equal to the hnal momentum, which is given by 
so the general energy momentum balance is given by 

d(E/j/c,pij) = -dm(7eC,7eu), (8) 

where u is the ejecta velocity in the center-of-momentum frame. The rocket’s 
total energy is given by Er = 'jRmRC^, while its total momentum is given by 
Ph = iRf^RV. Substituting these into Equation (|^ along with the signs for 
the 1-dimensional vectors and separating the equations gives the system 

dilRTHR) = -^edm, 
di'^RfURv) = '^eudm. 

The first of these equations represents the conservation of energy while the 
second represents the conservation of momentum. Forward m) at this point 
begins to evaluate the derivatives on the left hand side of the equation and 
upon substituting the fully relativistic form of u into the momentum equation 
arrives ultimately at Ackeret’s equation, unfortunately only by proceeding 
through a forest of differentiation and algebraic reduction. Suppose instead 
we substitute the hrst equation into the second to give 

di'yRmRv) = -ud{'yRmR). ( 10 ) 

The derivative on the left hand side of this equation is 

d{'yRmRv) = vd{jRmR) -F 'yRUiRdv, (11) 
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so defining 


/i = '-fRiriR, 


( 12 ) 


we may write this equation as 


{v + u)d^ + fidv = 0. 


(13) 


Equation (13) is a single, separable differential equation which describes any 


kind of rocket motion regardless of the velocity regime. Without specifying 
the expression to be used for u, we simply leave u = u{v) and perform the 
integration to hnd 

/i JO u' + u[v ) 

Since v{t = 0) = 0, we have 7 (t = 0) = 1 and thus no = ttiq. Therefore 


we may write equation (14), restoring the Lorentz factors and dehning the 


typical mass ratio (MR = minitiai/')^finai), as 


dv' 


= where VLiv) = T — , 

^ ' Jo v' + u{v') 


(15) 


Equation (15) forms the main result of the section and also the titular result 
of this work. In the next section we shall show how both the Tsiolkovsky 
equation and Ackeret’s relativistic equation emerge in the limit where u is 
treated Newtonianly and relativistically, respectively, as well as show that a 
series expansion of the integrand can be used to obtain relativistic corrections 
to a rocket which is flying in an intermediate regime. 


3. Evaluation of the Quadrature Solution in Varions Flight Regimes 


The evaluation of Equation (15) is accomplished by specifying a func¬ 
tional form for u{v) and performing the integration. We shall consider three 
such cases in this section corresponding in order to the Newtonian regime, 
relativistic regime, and the transition between the two. 


3.1. The Newtonian and Relativistic Regimes 

If V is far less than c, the relative velocity u is obtained through Newtonian 
addition of velocities, 

U = Ve — V, (16) 
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and the factor vjc is so small that 7 = 1 to any reasonable nnmber of 
signihcant hgnres. Thns we have 



dv' 

v' + {ve — v') 


V 

•) 

Ve 


(17) 


so 

MR = (18) 

which is Tsiolkovsky’s eqnation. 

Conversely, if v is close to c we must use the fully relativistic velocity 
addition formula, which is given by 


u = z ——where jSi = Vi/c. 




(19) 


Substituting this into Equation (|15|) and simplifying the expression gives 
1 


= V.l ) = a io - i ■ (20) 


dP' 


rP 13' 


P''- 


The rightmost integration is simple to perform and gives 


( 21 ) 


serving to effectively cancel out the factor of 7 in front of the exponential. 
The second integral is more difficult, but is found to be 

f3 d/3' _ 1 1 + /3 

Jo l-/3'2 “ 2 

which, when exponentiated, gives the familiar equation of Ackeret, 



MR 


/i+/5y^= 

0 -/ 3 / 


(23) 


This shows Equation (15) gives the correct results for both the relativistic 
and Newtonian limits. 
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3.2. The Transition to Relativistic Flight: A Series Solution 


Another result can be obtained from Equation (15) if we expand the 


integrand of Equation (22) in series, viz. 


1-/9^ 


= E/^ 


2u 


(24) 


u=0 


Integration of this series is trivial, producing 


/o 1-/5' 


1 .R oo oo q2i/+1 

d/3'= / d/3'E= E ^ 


72 


u=0 


„Vo2^ + l’ 


(26) 


which gives the expression for the relativistic mass ratio in terms of a product, 

1 /32^+^ \ 


oo / 1 o2p+1 ' 




n 


17=1 


/3e 2z/ + 1 


(26) 


If the product is taken to inhnity, the hnal result will be the same as Ack- 
eret’s equation. Indeed, it can be shown that if the leftmost product is 
reconstructed as the exponentiation of a series, this series converges exactly 


to the result of Equation (22), a trivial result given the simple steps taken 


to derive this form of the equation. 

Let us consider the case where /3 is large enough so that at least the hrst 
term of the rightmost product would not be unity, but not so large that an 
inhnite number of terms need be accounted for to successfully approximate 
the mass ratio. In this case it is appropriate to truncate the product in 


Equation (26) to a few low-order terms. Such truncations of this product 


show how the mass ratio varies as the rocket begins to enter the relativistic 
regime. Let us consider the family of rocket equations 


MR„ = e''/*V„(/?;A). 


(27) 


where 


n 

= n 

u=l 


/ 1 /32'^+i \ 

\K‘2iy + l) 


and po 


1 . 


(28) 


For the convenience of the reader, the hrst few Pn(/3; /3e) are given in Table 
1 in terms of n, Ug, and c and a plot of Pn(/3; /3e = 0.1) for n = 0,1, 2 and 5 
is shown in Figure 2. 

The set of equations denoted by MR,^ include the Tsiolkovsky equation 













Figure 2: A plot of Pn(/9; (^e = 0-1) for n = 0 (non-relativistic), 1, 2 and 5. 


(MRq) and Ackeret’s equation (MRoo). As the rocket transitions from the 
completely Newtonian regime to a weakly relativistic regime, the mass ratio 
stops being represented by MRq and begins being represented by MRi, given 
by 

MRi = /3e) = (29) 

The effect of this term is to require more propellant than would be required 
by the Tsiolkovsky equation alone to achieve the same hnal velocity were 
relativistic effects neglected. As the rocket moves still faster the second term 



Table 1: Relativistic corrections for the Tsiolkovsky equation up to 5th order. 
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in the product becomes significant and the mass ratio equation becomes 




(30) 


and so on until so many terms are required that it is simpler to use Ack- 



Figure 3: Mass Ratios versus f3 for various MR„ equations and j3e- The 
actual value of (3e does not change the relative morphology of the curves. 


eret’s equation. An obvious course of action from here is to determine the 
velocity range over which these low order truncations of the infinite product 
accurately predict the mass ratio. The difference between using MRi and 
Equation (23) for (3^ = 0.5 and [3 = 0.4 is about 0.5%, and as little as 0.3% 
for /3e = 0.999 and [3 = 0.4. At much larger values of /?, the first few terms 
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of this product are no longer sufficient, and it would be much easier to use 
the full relativistic equation. Figure 3 shows plots of In MR versus (3 given 
a hxed Vg for various n as well as the Tsiolkovsky equation with Ackeret’s 
equation as a limiting case. Note that the relative morphology of the curves, 
and thus the extent to which MR„ for a given n approximates the Ackeret 
equation, is invariant with respect to (3^. 


4 . Related Questions 

In the previous sections we showed the rocket equation can be reduced 
to a unique quadrature which is not dependent on the flight regime which is 
considered, and have used this fact to show explicitly the transition between 
Newtonian and relativistic flight regimes. In this section we discuss further 
questions which can be answered by the results we have developed, in par¬ 
ticular the case where Vg is relativistic but v is not and the case where the 
ejecta stream is massless but still has momentum. 

4.1. Relativistic Ejecta Stream on a Newtonian Rocket 

If f3g is considerable but /3 -C 1 we may write the equations of (|^ as 

dmR = —^gdm, 
d{mRv) = 'jgudm. 

This gives the differential equation 

udniR -|- d{mRv) = 0, 

which is indistinguisable from the classical case. In this case the gain in 
rocket momentum from the relativistic divergence of the ejecta momentum 
is exactly counteracted by the divergence in relativistic kinetic energy which 
the ejecta stream must gain before exiting the vehicle. This also explains 
the invariance of the relative morphology between the Ackeret equation and 
various MR„ with respect to jdg which was shown in the last section. 


(31) 

(32) 


4-2. Rocket Equation for A Massless Ejecta Stream 

Although the practicality of using a massless ejecta stream is disputed 
( 0 ), such cases have been considered in the literature before ( 0 . 1121 , 1111 ). 
A rocket equation governing this case may be derived from Equations ([^ and 
Equation (15) alike. For the analysis of this case we shall consider massive 
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reactants which are reacted in such a way as to create massless products that 
form the ejecta stream. This in particular ensures that the interpretation of 
the mass ratio remains obvious as the ratio of the total mass before and after 
the burn. 

In relativity the total energy of an object (rest-energy and kinetic energy) 
is given by 

E = yj (p ■ -1- {mc^y. (33) 

Basic quantum mechanical principles allow a second relationship between 
particle momentum and quantum frequency, 

E = hw, (34) 


so for a massless particle of quantum frequency u, 

hu 

p = —• 


(35) 


The existence of nontrivial energy and momentum relations implies it is 
possible in theory to create thrust using massless propellants alone, and thus 
there should be a well-dehned rocket equation governing such propellants. 
Furthermore it is known that massless particles (e.g. photons) travel at the 
speed of light regardless of the frame in which they are observed, and thus 


for the case of such particles u{v) = Ve = c. Equation (15) then gives 


rv Aqj' 

n{v)= / — = ln(l + /?), 
JO n' -f c 

which gives the rocket equation 


(36) 


MR = 


/1 + /3 

1-/3' 


(37) 


This same result may also be obtained by using the energy and momen¬ 


tum relations of Equations (34) and (35) and writing a system of equations 
similar to the equations 0 for a large number of massless particles. It 
is however more straightforward to perform the calculation using Equation 
(15). Furthermore this exactly agrees with Ackeret’s equation in the case 


where /3e = 1, an interesting result considering Ackeret’s equation is derived 
under the assumption of the ejecta stream having a hnite, nonzero mass. 
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5. Conclusions 

The dynamics involved in rocket motion are simple, and rely on basic 
conservation principles which are obeyed in every possible flight regime, in¬ 
cluding those which are not practically available due to technological limita¬ 
tions. This idea itself suggests that there should be a generic form the rocket 
equation, whose Newtonian and relativistic limits reproduce known results. 
We have shown that such a generic equation exists and can be simply de¬ 
rived from basic considerations. We have also used this generic equation to 
derive a series expansion of quasi post-Newtonian corrections which show 
the transition from non-relativistic to relativistic flight, and unexpectedly 
found that the traditional relativistic rocket equation is well approximated 
up to 0.5c by even the hrst-order relativistic correction to the Tsiolkovsky 
equation. Furthermore, we have found the difference in morphology between 
the traditional relativistic equation, Tsiolkovsky equation, and Tsiolkovsky 
equation with some number of relativistic expansion terms included is invari¬ 
ant with respect to the regime of the ejecta velocity. The most striking case 
of this is the hnding that a rocket traveling non-relativistically but utilizing 
a relativistic ejecta stream will not hnd any relativistic discrepancy, even in 
theory, between its true mass ratio and that predicted by the Tsiolkovsky 
equation. We have also found that the traditional relativistic equation will 
still accurately govern the mass ratio under the assumption that the ejecta 
stream is massless. 

It is possible that still further scenarios can be easily analyzed by for¬ 
mulating design-specihc expressions for u{v) and using the generic rocket 
equation we have derived to understand the mass ratio equation which re¬ 
sults from a given design. Examples of such scenarios might be hypothetical 
rockets using multiple ejecta streams of differing velocities, and perhaps op¬ 
erating within different velocity regimes. 
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